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Abstract. It is becoming clear that 'self-pollution' by the ejecta of massive asymptotic giant branch stars has 
an important role in the early chemical evolution of globular cluster stars, producing CNO abundance spreads 
which are observed also at the surface of unevolved stars. Observing that the ejecta which are CNO processed 
must also be helium enriched, we have modelled stellar evolution of globular cluster stars by taking into account 
this possible helium enhancement with respect to the primordial value. We show that the differences between the 
main evolutionary phases (main sequence, turn-off and red giants) are small enough that it would be very difficult 
to detect them observationally. However, the difference in the evolving mass may play a role in the morphology 
of the horizontal branch, and in particular in the formation of blue tails, in those globular clusters which show 
strong CNO abundance variations, such as M13 and NGC 6752. 
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J> 1. Introduction ter members of material nuclearly processed in their inte- 

rior, or pollution - partial or total - by external material 

^ ; The chemical inhomogeneities (spread in the abundances (^^^^ ^ ^gg^^ Cottrell and Da Costa 1981, 

^ . of CNO, O - Na and Mg - Al anticorrelation) observed m o'Antona et al. 1983, Kraft 1994, Cannon et al. 1998). The 

globular cluster (GC) members from the mam sequence confirmation that turn-off stars show peculiarities quite 

to the red giant branch impmge on problems such as clus- ^-^^^^ ^^^^^ observed in giants casts severe doubts on 

ter formation and early evolution (see, e.g., Kraft 1994, ^^^^ ^^^^ scenario (Gratton et al. 2001). 
Gratton et al. 2001). In some clusters the peculiarities are 

particularly wide spread and strong, as in M13 and NGC ^he primordial pollution hypothesis requires either 
6752 (Da Costa & Cottrell 1980, Norris et al. 1981, Smith contamination of intracluster material, out of which new 
& Norris 1993, Pilachowski et al. 1996, Kraft et al. 1997, ^tars form with a chemical composition different from a 
Sneden 2000, Gratton et al. 2001). The fact that these Preceding star generation, or partial contamination (su- 
clusters present also the "anomaly" of an extended and Perficial or deep) of pre-existing main sequence structures, 
exclusively blue horizontal branch has been noticed by The contaminating objects are generally identified with 
Catelan & de Freitas Pacheco (1995), who observe that intermediate mass asymptotic giant branch (AGB) stars 
clusters with super oxygen-poor giants have all very blue (4 " 6 Mq stars, see DAntona et al. 1983, Ventura et 
HBs. On the other hand, it seems most likely that there ^1. 2001, 2002; Thoul et al. 2002), which evolve rapidly 
is not a direct relation between the two features (i.e., blue enough to shed their envelopes in about 10^ yr. 
tail and anomalies), but rather that they may share the The detailed composition of these envelopes is not easy 
same origin, at least in part. to predict, as it depends on the uncertain physics of con- 
The most common hypotheses for the origin of the vection and mass loss in the intermediate mass AGB stars, 
chemical peculiarities are either deep mixing in the clus- but it will surely show an enrichment in helium, in part as 
a consequence of the second dredge-up, and in part due 



Send offprint requests to: V. Caloi to the third dredge-up: the helium mass fraction should 
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Fig. 1. Comparison of two isochrones of age 14 Gyr, Z=2 x 
10~^ and different helium abundance. The ZAHBs for the 
two Y values are also shown, down to a minimum HB mass 
of M=0.6 Mq. 

be ^ 0.29 instead of ^^0.24, the cosmic value (Ventura 
et al. 2001)0. When these envelopes are expelled at low 
velocity by the stellar bodies during the AGB evolution, 
various consequences can arise. Given the depth of the 
potential well, a significant fraction of this material is not 
lost by a globular cluster, unless there is a significant in- 
teraction with external systems (e.g. another cluster, or a 
giant molecular cloud, or the galactic bulge). New stars 
can form out of this matter, or pre-existing diffuse mate- 
rial can be polluted, or low mass unevolved stars can be 
polluted to various degrees. 

In order to evaluate the impact of a possible popula- 
tion having a different chemical composition on the color 
magnitude diagram of a globular cluster, in this article 
we compare isochrones having primordial helium abun- 
dance 0, which we set at Y=0.23 for Z=2xl0^"', and at 
Y=0.24 for Z^IO"'^, with models having a main sequence 

^ Notice that the value Y=0.29 for the ejecta of the most 
massive AGBs is substantially conservative, as the Ventura et 
al. models make very conservative assumptions on the modal- 
ities of the third dredge up. 

We take Y=0.24 as standard value for the metallicity 
Z=10~^. Implicitly, we are assuming that the Big Bang value 
was Y=0.23, and that there is a small increase in Y for the 
matter reaching a metallicity Z=10~'^. Alternatively, we may 
regard Y=0.23 and Y=0.24 both as values indicative of the 
helium abundance emerging from the Big Bang. The value 
Y=0.228±0.005 was derived by Pagel et al. (1992), whereas 
today a larger value, Y=0.244±0.002 results from Izotov and 
Thuan (1998) analysis of the helium emission lines in HII re- 
gions of metal poor dwarf galaxies. 



Fig. 2. The same comparison as in Fig. for Z=10 ^ . In 
this case, we have extended the ZAHB down to masses 
as small as M=0.5 Mofor Y=0.24 and M=0.491 M© for 
Y=0.28. These additional (bluer) models below 0.6 Mq 
are necessary for the simulations presented in Sect. 2.1. 

enhanced helium Y=0.28. We compute models for all evo- 
lutionary phases and compare the location of the main 
sequence (MS), turn-off (TO), red giants (RGB) and hor- 
izontal branch (HB). The lower main sequence models, 
and the detailed evolution along the RGB are considered 
in a forthcoming paper (Montalban et al. 2002). In Figs. 
1 and 2, the low mass cutoff of each main sequence is 0.6 
Mq. Notice that the Y=0.28, 0.6 Mq models are much 
more luminous (by more than 0.4 mag) than the standard 
helium models. Nevertheless, the relative locations of the 
standard and enhanced helium isochrones show that it 
is not obvious how to discriminate between a population 
with uniform helium content and one with a spread in he- 
lium, for what concerns TO and RGB evolution. On the 
other hand, the morphology of the HB is affected by a he- 
lium variation, directly through the difference in the TO 
mass for a given age and we argue that a helium variation 
may be the main reason for the occurrence of extreme blue 
tails in clusters such as NGC 6752 and M13. 

This conclusion differs from the common hypothesis 
put forward in recent years: a helium enrichment in the 
envelope as a consequence of deep mixing in today evolv- 
ing red giants has in fact been suggested as a way to ob- 
tain blue and very blue HBs (Langer & Hoffman 1995, 
Sweigart 1997). As a matter of fact, an increased helium 
abundance can help in reaching a HB position bluer than 
the RR Lyrae variables (at the expense of a more or less 
substantial increase in luminosity), but the extremely blue 
locations always require an extreme mass loss, because the 
envelope mass has to be in any case smarter than 0.01 Mq 
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(see the discussion in Caloi 2001). In addition, Weiss et al. 
(2000) have shown that the hehum enhancement needed 
to explain the exclusively blue HB morphology would pro- 
duce related chemical inhomogeneities in red giants much 
stronger than observed. So an envelope helium content 
larger than the cosmological one does not help in obtan- 
ing the structures which make up the long blue tails in, 
e.g., M13 and NGC 6752. We show that the situation is 
different if helium is larger in the whole stellar structure, 
that is, if cluster stars form out of material with enhanced 
helium. We considered this hypothesis at the light of the 
present state of the art in stellar modelling. 



2. The models and the comparison with 
observations 

Appropriate evolutionary models have been computed 
with the code ATON2.0, described in Ventura et al. (1998) 
and Mazzitelli et al. (1999). We considered models with 
Y=0.28, and two values for Z: 2 lO"* and 10"^, to be com- 
pared with tracks with Y=0.23, Z=2 10""*, and Y=0.24, 
Z=10~^. The increase in the helium content has been 
taken meaningful, but not such as to a priori alter dra- 
matically the currently accepted evolution of a GC low 
mass star, and it roughly agrees with the expected he- 
lium enhancement in the ejecta of massive, metal poor 
AGB stars. We have built isochrones for GC ages, and 
transformed the luminosity and Toff's into observational 
magnitudes and colors by using the transformations by 
Castelh et al. (1997) .| 

Figures |l| and || show the location of the 14 Gyr 
isochrones in the HR diagram, and the relevant varia- 
tions in the main loci. Similar variations hold for differ- 
ent ages. The computed isochrones show that there is a 
small age difference for the same TO magnitude ('--^0.8 
Gyr for Z=10~^, and ~ 1 Gyr for Z=2xl0~'* respectively, 
the larger helium the younger), while for the same age 
one finds 0.05 (0.07) mag difference (the larger helium the 
fainter, see Figs. and |). The main sequence location for 
the higher helium is 0.12 (0.13) mag fainter and the he- 
lium core mass at the red giant tip is 0.008 smaller for 
the higher helium tracks (Mcorc = 0.497 and 0.489 Mq, re- 
spectively, for Z=10^'^). The color difference between the 
low and high helium isochrones is ^ 0.015 mag both along 
the RGB and ~1 mag below the MS TO. In a cluster with 
a dispersion in helium abundance, such a small difference 
is certainly beyond present detection capabilities. A small 
population of binaries would make even more difficult to 
distinguish a possible broadening of the mean loci. 



^ Tables of isochrones and HB models in Mv, B-V and V-I 
are available on request. We are working to make available the 
models also in other (HST) passbands, as well as a computation 
tool to build HB synthetic populations, taking into account 
helium content spreads. 



Table 1. Input data for the HB models 
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2.1. The horizontal branch 

Table 1 shows the input parameters for the HB models, 
namely the helium core mass and the surface helium 
abundance Yq. Mc is the helium core mass at the RG tip, 
appropriate for ages between 10 and 16 Gyr, according to 
the present models. Yc is increased with respect to the 
initial value, due to the first dredge up. The zero age hor- 
izontal branch (ZAHB) loci are shown in Figs. |l| and |[ 
while Fig. ^ shows the HB mass versus the ZAHB B-V 
for the Z==10^^ case. 

For what concerns the actual distribution along the 
HB locus, the relevant feature is that the isochrone with 
initial helium content of Y=0.28 reaches the giant tip with 
a mass smaller of 0.05 Mq than in the case with standard 
helium. This is a well known result of the evolution — see, 
e.g., Iben and Renzini (1984) — which in our case provides 
an important clue to the HB morphology of those globular 
clusters in which we may suspect that there is a spread 
or a dichotomy in the helium abundance. Although Fig. |^ 
shows that, at a given B-V color, the Y=0.28 ZAHB mass 
is smaller by ~ O.OIMq than the mass for lower helium, 
we expect that the helium enriched stars will occupy bluer 
HB loci, mainly due to the difference in the evolving total 
mass. 

In order to make this point more clear, we computed 
a few simulations of the HB distributions expected for 
different assumptions on mass loss and helium content of 
the sample stars. These simulations are not intended to 
cover in an exhaustive way all possible cases, but sim- 
ply to exemplify the main new point of the present work. 
We make the assumption that the evolving giants include 
two components: the first one has initial abundance of 
Z=:0.001 and Y=0.24; the second one has the same metal 
abundance, but a uniform distribution of the initial he- 
lium content between Y=0.24 and Y=0.28. We also ran 
other simulations testing the hypothesis of a gap in the 
helium content distribution between the two components. 
This would be the case if a new generation of stars were 
born in the cluster directly from the wind ejecta of the 
massive AGBs. 

We assume that the individual HB masses are the re- 
sult of a mass loss process which goes on along the RGB, 
with a Reimers' (1997) like mass loss law. We assume that 
the average mass lost in this process is 5M, with a gaussian 
distribution having mass dispersion a (in this way, stars 
with identical chemical compositions may populate differ- 
ent part of the ZAHBs). Further, as Reimers' mass loss 
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Fig. 3. The mass versus B-V distribution in ZAHB for 
the two helium contents in the Z—IQ^^ case. Formally, the 
Y=0.28 models are rerfrfer that the Y=0.24 models, for the 
same total mass, except for the smaller masses. However, 
the mass loss mechanism at the red giant branch operates 
on a smaller total mass for the high helium case, so that 
it is easier to achieve a smaller remnant mass and bluer 
stars. 

rate increases with decreasing Mrg (see the discussion in 
Lee et al. 1994), we also take into account that more mass 
is lost on average from the helium rich (and smaller mass) 
evolving giants. We ran various simulations with different 
values of 6M and oQ, and computed the distribution of 
stars along the HB for several ages. Here we show a few 
examples of the results obtained for an age of 13 Gyr. Very 
similar results can be obtained for larger (smaller) ages, 
by decreasing (increasing) accordingly the average mass 
loss. The total number of HB stars originated from main 
sequence stars with Y=0.24 is taken to be 250. To these, 
we add from 100 to 150 stars deriving from main sequence 
structures with helium content varying from 0.24 to 0.28, 
as mentioned before. 

In the four panels of Fig. ^, the HB members with origi- 
nal Y=0.24 are indicated by open triangles, while squares 
indicate stars deriving from the He-rich population de- 
scribed above (original Y from 0.24 to 0.28); evolution 
out of the ZAHB is taken into account, together with a 
small spread in V (ctv — 0.005) and in B-V (ctb-v = 0.01), 
which again simulate the impact of observational errors. 
Cluster members with low helium abundance give origin 
to the bulk of HB distribution, while, as expected, the 

We do not assume any explicit dependence of mass loss on 
the metallicity, as we are only giving examples of the resulting 
distributions for one metallicity only. 




-0.4 -0.3 -0.2 -0.1 0.1 

B-V 

Fig. 5. A simulation is shown together with the corre- 
sponding evolutionary tracks. 



high helium structures are able to populate much bluer 
HB regions. This is in part a consequence of the smaller 
mass evolving at the TO (up to 0.05 M© for the stars with 
Y=0.28), and in part because of the extra mass loss due 
to the smaller mass of the evolving giants (up to 0.03 Mq 
for Y=0.28; see Lee et al. 1994). 

When the average mass loss is small, most stars would 
evolve on the red side of the HB (top panels of Fig. 
However, stars on the red HB and within the RR Lyrae 
strip computed with a distribution of initial values of 
the helium abundance would be on average more lumi- 
nous and yield a substantially thicker distribution than 
expected in the standard (low helium) case. For clusters 
more metal rich than considered here, in which the HB 
is restricted to the red of the RR Lyrae gap, this fea- 
ture has probably been observed: in 47 Tuc (Briley 1997), 
CN-strong HB stars appear on the average slightly more 
luminous (0.05 mag) and tend to be bluer, than their CN- 
weak counterparts. By increasing the average mass loss in 
our simulations, (top right panel in Fig. ^ the helium rich 
stars begin to populate the whole HB extension from blue 
to red, even assuming a relatively small dispersion in the 
mass loss of cr = 0.012. 

If the HB is of the type shown in the bottom left panel 
of Fig. the luminosity level of the HB is determined by 
the Y=0.24 stars, so that the problem of predicting too 
high a luminosity for the HB if there is a helium spread is 
not relevant: all the helium rich stars, in fact, are confined 
to the blue side of the HB, due to their smaller mass and 
larger mass loss in the RGB phase. 

In exclusively blue HBs (bottom right of Fig. ||), again 
the higher helium structures are confined on the blue tail. 
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Fig. 4. The four panels show the synthetic HB for different choices of the average mass loss (dM in Mq), and dispersion 
in mass loss (ct, in M0). The age assumed is 13 Gyr and the metallicity is Z=10~^. The simulations contain 400 stars, 
of which 250 have Y=0.24 (triangles) and the remaining 150 have Y randomly chosen between 0.24 and 0.28 (squares). 



where helium abundance does not affect noticeably ZAHB 
position and early evolution (see, e.g., Caloi 2001). 

In Fig. ^ we show the typical HB distribution in a 
simulation including 250 stars with Y=0.24 and 150 stars 
with Y uniformely distributed between 0.24 and 0.28. The 
corresponding tracks are also shown. The chosen mass loss 
rate is such that only the helium rich stars populate the 
extreme blue tail. Other experiments with different values 
of the mean mass loss and mass dispersion, are shown in 
Figs. H and 0. Here we have also introduced a gap in the 
helium range of the enriched population (Y from 0.26 to 
0.28 in Fig. | and a fixed Y=0.28 for Fig. ^. In Figs. 
^ and 1^ we also show the histogram of the luminosity 
function of HB stars. We see that a dip in the distribution 
appears at the magnitude corresponding to the change 
in the helium abundance, and we consider this to be a 
meaningful feature. Other features at larger Mv do not 
appear to be statistically significant. 



With values for the parameters largely in the accepted 
range, we can obtain the morphologies of well known 
globular clusters with blue tails (cfr. Piotto et al. 1999, 
Buonanno et al. 1986, Ferraro et al. 1998). A very blue 
clump well separated from a "normal" blue HB (that is, 
derived from giants with standard helium) can be obtained 
if helium in the enriched population peaks at Y=0.28 (Fig. 
^). The bluest HB structures have a hydrogen rich enve- 
lope of the order of 0.001 M0, and evolve vertically in 
the CM diagram at a constant rate (as predicted by all 
evolutionary computations), partially filling the gap. 



3. Conclusion 

The problem of chemical inhomogeneities in GCs is part 
of the wider question of their chemical history, which ap- 
pears more and more complex, and differing from clus- 
ter to cluster. As put in evidence, e.g., by Sneden (2000), 
the neutron-capture elements Ba and Eu have variable 
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Fig. 6. Simulation containing 380 stars, 250 of which have 
Y=0.24, and 130 hehum randomly chosen between 0.26 
and 0.28. 
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Fig. 7. Simulation containing 350 stars, 250 of which have 
Y=0.24, and 100 have Y=0.28. 

abundances without obvious connection to overall cluster 
metallicity. Besides, variations in Si content are also ob- 
served from cluster to cluster. 

The fact that field stars do not show evidence of 
the substantial abundance re-shuffling between CNO el- 
ements and Na, Mg and Al strongly suggests some intra 
cluster processes at the origin of the phenomenon. Here we 



examine one possible aspect of these (admittedly) rather 
vague processes. 

We discussed the evolutionary properties of GC stars 
in which variations of the primordial helium content are 
allowed. These variations are attributed to self-pollution 
by an initial population of intermediate mass stars, and 
are expected on the basis of the same stellar models which 
are thought to produce the CNO variations among mem- 
bers of several GCs. 

The observation of chemical peculiarities (dichotomy 
in the strength of CN lines, Na and Al enhancement) 
in some GC giants suggested already a long time ago 
(CottreU and Da Costa 1981, Norris et al. 1981) that in- 
termediate mass stars may have played a role in creat- 
ing some inhomogeneity in the primordial cloud of certain 
GCs. In particular, Norris et al. (1981) put in relation 
chemical peculiarities with blue tails (NGC 6752), and 
suggested precisely the mechanism which is at the basis 
of the present investigation. 

We find that a spread in the helium content does not 
affect the morphology of the MS, TO and RGB in an eas- 
ily observable way. The simulations in Figs. § - show 
that also the HB luminosity level is affected in a limited 
way, since only an RR Lyrae region populated exclusively 
by structures with the maximum Y allowed (0.28) would 
stand out clearly as peculiarly luminous. On the other 
hand, such an occurrence should not be verified with stan- 
dard values for the mean mass loss (AM ^ 0.2 M©, see 
Lee et al. 1994). 

So age determination should not be affected in a sub- 
stantial way. There may be a small effect if we use as 
distance indicator the subdwarf main sequence, since it is 
reasonable to assume that subdwarfs have a helium con- 
tent not affected by self-pollution. 

On the other hand, the helium spread may constitute 
part of the "second parameter" problem. The presence of 
very blue extensions in many clusters with predominantly 
blue HB, and sometimes also in HBs with RR Lyrae vari- 
ables and red members -NGC 2808 (Walker 1999), M62 
(Caloi et al. 1987, Brocato et al. 1996)- may be related to 
the second parameter problem, but this side of the ques- 
tion will not be considered here. The extremely small enve- 
lope masses necessary to populate these blue tails require 
mass losses, and/or ages, much larger than the average 
values required to explain also exclusively, but not ex- 
tended, blue HBs, in case of uniform cosmological helium 
abundance. As we saw before, the smaller TO masses (for 
a given age) of a helium enriched structure are crucial to 
obtain the tiny envelopes of very blue HB stars. 

The difficulties with the present hypothesis are mainly 
due to the enrichment scenario, since it is not easy to pol- 
lute a large fraction of original intraclustcr matter in a 
substantial way (Cannon et al. 1998). We note that the 
requirement of the pollution of a large fraction of the clus- 
ter population comes from the roughly bimodal distribu- 
tion of CN abundances and from the well developed blue 
tails in many clusters, if they are related to the helium 
enhancement phenomenon (M13, NGC 6752). 
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In the self-enrichment scenarios, if the extreme blue 
tails can indeed be attributed to the difference among 
the total masses of the stars evolving in the RGB, due 
to their different helium content, the helium content must 
be larger also in the center of the stars, to affect their 
hydrogen burning lifetimes, and therefore we should con- 
clude that the self-pollution has occurred cither on the gas 
from which the helium rich stellar generation was formed, 
or it has occurred during the phases of the stars evolution 
in which they were fully convective. It is also possible that 
the ejecta of AGBs, collected at the center of the cluster, 
directly form other stellar generations. If this is the case, 
some gaps in the HB stellar distributions could rc^flect dis- 
continuities in the helium content of the constituent mat- 
ter, due to the difference in the helium abundance in the 
ejecta of AGBs of different mass (Ventura (it al. 2002) and 
to some peculiar modalities of star formation from these 
ejecta. 
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